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A description of the intermediate and high energy hadronic interaction models used in the FLUKA 
code is given. Benchmarking against experimental data is also reported in order to validate the 
model performances. Finally the most recent developments and perspectives for nucleus-nucleus 
interactions are described together with some comparisons with experimental data. 



I. GENERALITIES 

FLUKA 0, H H A is a multipurpose trans- 
port Monte Carlo code, able to treat hadron-hadron, 
hadron-nucleus, neutrino, electromagnetic, and /j, in- 
teractions up to 10000 TeV. Charged particle trans- 
port (handled in magnetic field too) includes all rele- 
vant processes 1]. About nucleus-nucleus collisions, 
since ion-ion nuclear interactions were not yet treated 
in FLUKA, past results have been obtained in the su- 
perposition model approximation, where primary nu- 
clei (0-10000 TeV/A) were split into nucleons before 
interacting. With the integration of ion interaction 
codes (see Section ITVj) and the cross section parame- 
terization, this approximation is now obsolete. 

FLUKA is based, as far as possible, on original and 
well tested microscopic models. Due to this "micro- 
scopic" approach to hadronic interaction modelling, 
each step is self-consistent and has solid physical 
bases. Performances are optimized comparing with 
particle production data at single interaction level. No 
tuning whatsoever is performed on "integral" data, 
such as calorimeter resolutions, thick target yields, 
etc. Therefore, final predictions are obtained with a 
minimal set of free parameters, fixed for all energies 
and target/projectile combinations. Results in com- 



plex cases as well as scaling laws and properties come 
forth naturally from the underlying physical models 
and the basic conservation laws are fulfilled a priori. 



II. HADRON-NUCLEON INTERACTION 
MODELS 

A comprehensive understanding of hadron-nucleon 
(h-N) interactions over a wide energy range is of 
course a basic ingredient for a sound description of 
hadron-nucleus ones. Elastic, charge exchange and 
strangeness exchange reactions are described as far as 
possible by phase-shift analysis and/or fits of exper- 
imental differential data. Standard eikonal approxi- 
mations are often used at high energies. 

At the low energy end (below 100 MeV) the p-p 
and p-n cross sections are rapidly increasing with de- 
creasing energy. The n-p and the p-p cross sections 
differ by about a factor three at the lowest energies, as 
expected on the basis of symmetry and isospin consid- 
erations, while at high energies they tend to be equal. 

The total cross section for the two isospin com- 
ponents present in the nucleon-nucleon amplitude is 
given by: 



O-Q 



— vpp 

2(7 n p 0~pp . 



*on leave from GSI, Darmstadt, Germany (Present address: 
SLAC, P.O. Box 4349, MS 48, Stanford, CA 94309, USA) 
ton leave from INFN, Sezione di Milano, via Celoria 16, 1-20133, 
Milan, Italy 



The same decomposition can be shown to apply for 
the elastic and the reaction cross sections too. 

The cross section for pion-nucleon scattering is 
dominated by the existence of several direct reso- 
nances, the most prominent one being the A(1232). 
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Given the pion isotopic spin (T = 1), the three 7r 
charge states correspond to the three values of T z . 
Thus, in the pion-nucleon system two values of T are 
allowed : T = ^ and T = §, and two independent 

scattering amplitudes, Ai and A3 , enter in the cross 

° 2 2 

sections. Using Clebsch-Gordan coefficients all dif- 
ferential cross sections can be derived from the three 
measured ones: a (i: + p — > 7r + p), a {n~p — » 7r~p), and 
a VK~p — > 7r°n) . 

As soon as particle production is concerned (inelas- 
tic hadron-nucleon interactions), the description be- 
comes immediately more complex. Two families of 
models are adopted, depending on the projectile en- 
ergy: those based on individual resonance production 
and decays, which cover the energy range up to 3- 
5 GeV, and those based on quark/parton string mod- 
els, which can provide reliable results up to several 
tens of TeV. 



B. Inelastic h-N interactions at high energies 
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A. h-N interactions at intermediate energies 

The inelastic channel with the lowest threshold, sin- 
gle pion production, opens already around 290 MeV in 
nucleoli- nucleon interactions, and becomes important 
above 700 MeV. In pion-nucleon interactions the pro- 
duction threshold is as low as 170 MeV. Both reactions 
are normally described in the framework of the isobar 
model: all reactions proceed through an intermedi- 
ate state containing at least one resonance. There 
are two main classes of reactions, those which form 
a resonant intermediate state (possible in 7r-nuclcon 
reactions) and those which contain two particles in 
the intermediate state. The former exhibit bumps in 
the cross sections corresponding to the energy of the 
formed resonance. Examples are reported below: 

Ni + N 2 -> N[ + A(1232) N{ + + tt 
tt + N -> A(1600) -► tt' + A(1232) -> 

-> 7r'+TT"+N' 

Ni + N 2 -> Ai(1232) + A 2 (1232) -> 

N[ + Ttl + N' 2 + 7T 2 

Partial cross sections can be obtained from one- 
boson exchange theories and/or folding of Breit- 
Wigner with matrix elements fixed by N-N scatter- 
ing or experimental data. Resonance energies, widths, 
cross sections, and branching ratios are extracted from 
data and conservation laws, whenever possible, mak- 
ing explicit use of spin and isospin relations. They 
can be also inferred from inclusive cross sections when 
needed. 

For a discussion of resonance production, see for 
example HIE0- 



FIG. 1: Leading two-chain diagram in DPM for p — p scat- 
tering. The colour (r^red, b— >blue, g^green, r — >antired, 
b ^antiblue and g — >antigreen) and quark combination 
shown in the figure is just one of the allowed possibilities. 
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FIG. 2: One of leading two-chain diagrams in DPM for 
7r + — p scattering. The colour and quark combination 
shown in the figure is just one of the allowed ones. 

As soon as the projectile energy exceeds a few GeV, 
the description in terms of resonance prodution and 
decay becomes more and more difficult. The number 
of resonances which should be considered grows expo- 
nentially and their properties are often poorly known. 
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pi+ p -> CVC- X, 250 GeV/c 
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X F _ 

FIG. 3: Feynman x* F spectra of positive and negative 
particles from (ir + , p) at 250 GeV/c . Exp. data (symbols) 
from 0. 



pi+ + p -> C+/C- X, 250 GeV/c 




FIG. 4: Transverse momentum (p t ) spectra of positive and 
negative particles from (ir + , p) at 250 GeV/c . Exp. data 
(symbols) from |g. 



Furthermore, the assumption of one or two resonance 
creation is unable to reproduce the experimental find- 
ing that most of the particle production at high ener- 
gies occurs neither in the projectile nor in the target 
fragmentation region, but rather in the central region, 
for small values of Feynman x variable. Different mod- 
els, based directly on quark degrees of freedom, must 
be introduced. 

The features of "soft" interactions (low-px interac- 
tions) cannot be derived from the QCD Lagrangian, 
because the large value taken by the running cou- 
pling constant prevents the use of perturbation the- 
ory. Models based on interacting strings have emerged 
as a powerful tool in understanding QCD at the soft 
hadronic scale, that is in the non-perturbative regime. 
An interacting string theory naturally leads to a topo- 
logical expansion. The Dual Parton Model (DPM) 
is one of these models and it is built introducing par- 
tonic ideas into a topological expansion which explic- 
itly incorporates the constraints of duality and uni- 
tarity, typical of Regge's theory. In DPM hadrons are 
considered as open strings with quarks, antiquarks or 
diquarks sitting at the ends; mesons (colourless com- 
bination of a quark and an antiquark qq) are strings 
with their valence quark and antiquark at the ends. 
At sufficiently high energies the leading term in the in- 
teractions corresponds to a Pomeron (IP) exchange (a 
closed string exchange), which has a cylinder topology. 
When an unitarity cut is applied to the cylindrical 



Pomeron, two hadronic chains are left as the sources 
of particle production. As a consequence of colour ex- 
change in the interaction, each colliding hadron splits 
into two coloured system, one carrying colour charge c 
and the other c. The system with colour charge c (c) of 
one hadron combines with the system of complemen- 
tary colour of the other hadron, to form two colour 
neutral chains. These chains appear as two back-to- 
back jets in their own centre-of-mass systems. 

The exact way of building up these chains depends 
on the nature of the projectile-target combination 
(baryon-baryon, meson-baryon, antibaryon-baryon, 
meson- meson) : examples are shown in figs. Q and y] 
Further details can be found in the original DPM ref- 
erences or in 0| . 

The chains produced in an interaction are then 
hadronized. DPM gives no prescriptions on this stage 
of the reaction. All the available chain hadronization 
models, however, rely on the same basic assumptions, 
the most important one being chain universality, that 
is chain hadronization does not depend on the par- 
ticular process which originated the chain, and un- 
til the chain energy is much larger than the mass of 
the hadrons to be produced, the fragmentation func- 
tions (which describe the momentum fraction carried 
by each hadron) are the same. As a consequence, frag- 
mentation functions can in principle be derived from 
hard processes and e + e~ data and the same functions 
and (few) parameters should be valid for all reactions 
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and energies; actually mass and threshold effects are 
non-negligible at the typical chain energies involved 
in hadron-nucleus reactions. Transverse momentum 
is usually added according to uncertainty considera- 
tions. The examples in figs. |3| and 0] show the ability 
of the FLUKA model, based on DPM, to reproduce 
the features of particle production; further examples 
can be found in 0, . 



III. MAIN STEPS OF HADRON-NUCLEUS 
INTERACTIONS 

High energy hadron-nucleus (h-A) interactions can 
be schematically described as a sequence of the fol- 
lowing steps: 

• Glauber-Gribov cascade 

• (Generalized) IntraNuclear cascade ((G)INC) 

• Preequilibrium emission 

• Evaporation/Fragmentation/Fission and final 
deexcitation 



The Glauber-Gribov cascade and the 
formation zone 




FIG. 5: Leading two-chain diagrams in DPM for p — 
A Glauber scattering with 4 collisions. The colour and 
quark combinations shown in the figure are just one of the 
allowed possibilities. 



The Glauber formalism [T(l ITlj provides a powerful 
and elegant method to derive elastic, quasi-elastic and 
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FIG. 6: Multiplicity distribution of negative shower par- 
ticles for 250 GeV/c K + on aluminium and gold targets. 
Data from |Ts||. 



absorption h-A cross sections from the free h-N cross 
section and the nuclear ground state only. Inelastic 
interactions are equivalent to multiple interactions of 
the projectile with v target nucleons. The number of 
such "primary" interactions follows a binomial distri- 
bution (at a given impact parameter, b): 



P rv (b) 



A 



P r "(b) [l-P r (b)} 



A-v 



where P r (b) = UhNr Tr {b) , and T r (b) is the profile 
function (folding of nuclear density and scattering pro- 
files along the path). On average: 



O'hAabs(S) 



< V > = 



d 2 b 



Za h 



ip r 



No hr . 



0~hAabs 



1 - (1 - o- hNr (s)T r (b)Y 



The Glauber-Gribov model Il3 . repre- 

sents the diagram interpretation of the Glauber cas- 
cade. The v interactions of the projectile orig- 
inate 2v chains, out of which 2 chains (valence- 
valence chains) struck between the projectile and tar- 
get valence (di)quarks, 2{v — 1) chains (sea-valence 
chains) between projectile sea q — q and target valence 
(di)quarks. 

A pictorial example of the chain building process is 
depicted in fig. El for p-A: similar diagrams apply to 
7T-A and p -A respectively. 
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FIG. 7: Shower, grey, and black tracks multiplicities for 
7r~ on emulsion, as a function of the projectile momentum. 
Open symbols are experimental data from various sources, 
full symbols are FLUKA results. 



The distribution of the projectile energy among 
many chains naturally softens the energy distributions 
of reaction products and boosts the multiplicity with 
respect to hadron-hadron interactions. The building 
up of the multiplicity distribution from the multi- 
ple collisions can be appreciated from fig. |SJ where 
the multiplicity distribution for Al and Au targets at 
250 GeV/c are presented together. In this way, the 
model accounts for the major A-dependent features 
without any degree of freedom, except in the treat- 
ment of mass effects at low energies. 

The Fermi motion of the target nucleons must be in- 
cluded to obtain the correct kinematics, in particular 
the smearing of pr distributions. All nuclear effects on 
the secondaries are accounted for by the subsequent 
(G)INC. 

The formation zone concept is essential to under- 
stand the observed reduction of the re-interaction 
probability with respect of the naive free cross sec- 
tion assumption. It can be understood as a "materi- 
alization" time. At high energies, the "fast" (from the 
emulsion language) particles produced in the Glauber 
cascade have a high probability to materialize already 
outside the nucleus without triggering a secondary 
cascade. Only a small fraction of the projectile energy 
is thus left available for the INC and the evaporation. 

The Glauber cascade and the formation zone act 
together in reaching a regime where the "slow" part 

MOMT005 



of the interaction is almost independent of the particle 
energy. This can be easily verified looking at charged 
particle average multiplicities and multiplicity distri- 
butions as a function of energy (fig. 0). "Fast" (or 
"shower") tracks (charged particles with [3 = - > 
0.7), coming from the projectile primary interactions, 
show the typical « logarithmic increase observed for 
h-N interactions. As shown already in fig. El the av- 
erage multiplicity and its variance are directly related 
to the distribution of primary collisions as predicted 
by the Glauber approach. Due to the very slow varia- 
tion of h-N cross section from a few GeV up to a few 
TeV, the Glauber cascade is almost energy indepen- 
dent and the rise in the multiplicity of "fast" particles 
is related only to the increased multiplicity of the el- 
ementary h-N interactions. 

Due to the onset of formation zone effects, most of 
the hadrons produced in the primary collisions escape 
from the nucleus without further reinteractions. Fur- 
ther cascading only involves the slow fragments pro- 
duced in the target fragmentation region of each pri- 
mary interaction, and therefore it tends quickly to sat- 
urate with energy as the Glauber cascade reaches its 
asymptotic regime. This trend is well reflected in the 
average multiplicity (and multiplicity distribution) of 
"gray" tracks (charged particles with 0.3 < (3 < 0.7), 
which are mostly protons produced in secondary col- 
lisions during the INC and preequilibrium phases. 

At the end of the cascading process, the residual 
excitation energy is directly related to the number 
of primary and secondary collisions which have taken 
place. Each collision is indeed leaving a "hole" in the 
Fermi sea which carries an excitation energy related 
to its depth in the Fermi sea. Evaporation products, 
as well as residual excitation functions, should reach 
an almost constant condition as soon as the Glauber 
mechanism and the formation zone are fully devel- 
oped. This can actually be verified by looking at the 
production of "black" tracks (charged particles with 
[3 < 0.3), which are mostly evaporation products. The 
data reported in fig.[7|do indeed demonstrate how they 
saturate as well, and how this property is well repro- 
duced on the basis of the outlined assumptions. 

B. (Generalized) IntraNuclear cascade 

At energies high enough to consider coherent effects 
as corrections, a h-A reaction can be described as a 
cascade of two-body interactions, concerning the pro- 
jectile and the reaction products. This is the mecha- 
nism called IntraNuclear Cascade (INC). INC models 
were developed already at the infancy of the computer 
era with great success in describing the basic features 
of nuclear interactions in the 0.2-2 GeV range. Mod- 
ern INC models had to incorporate many more ideas 
and effects in order to describe in a reasonable way 
reactions at higher and lower energies. Despite parti- 
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cle trajectories are described classically, many quan- 
tistic effects have to be incorporated in these (G)INC 
models, like Pauli blocking, formation time, coherence 
length, nucleon antisymmetrization, hard core nucleon 
correlations. A thorough description of the (G)INC 
model used in FLUK A can be found in P, 0] ■ 

C. Preequilibrium 

At energies lower than the 7r production thresh- 
old a variety of preequilibrium models have been de- 
veloped [lf| following two leading approaches: the 
quantum-mechanical multistep model and the exci- 
ton model. The former has a very good theoretical 
background but is quite complex, while the latter re- 
lies on statistical assumptions, and it is simple and 
fast. Exciton-based models are often used in Monte 
Carlo codes to join the INC stage of the reaction to 
the equilibrium one. 

In the FLUKA implementation the INC step goes 
on until all nucleons are below a smooth threshold 
around 50 MeV, and all particles but nucleons (typ- 
ically pions) have been emitted or absorbed. At the 
end of the INC stage a few particles may have been 
emitted and the input configuration for the preequi- 
librium stage is characterized by the total number of 
protons and neutrons, by the number of particle-like 
excitons (nucleons excited above the Fermi level) , and 
of hole-like excitons (holes created in the Fermi sea by 
the INC interactions), and by the nuclear excitation 
energy and momentum. All the above quantities can 
be derived by properly recording what occurred dur- 
ing the INC stage. The exciton formalism of FLUKA 
follows that of M. Blann and coworkers [l7l TTM fl9l . l20 | . 
with some modifications detailed in H. 



D. Evaporation, fission and nuclear break-up 

At the end of the reaction chain, the nucleus is a 
thermally equilibrated system, characterized by its ex- 
citation energy. This system can "evaporate" nucle- 
ons, fragments, or 7 rays, or can even fission, to dis- 
sipate the residual excitation. 

Neutron emission is favoured over charged particle 
emission, due to the Coulomb barrier, expecially for 
medium-heavy nuclei. Moreover, the excitation en- 
ergy is higher in heavier nuclei due to the larger cas- 
cading chances and to the larger number of primary 
collisions in the Glauber cascade at high energies. The 
level density parameter a ~ A/8 MeV _1 is higher too, 
thus the average neutron energy is smaller. Therefore, 
the neutron multiplicity is higher for heavy nuclei than 
for light ones. 

The FLUKA evaporation module is based on the 
standard Weisskopf-Ewing formalism 21]. Latest 
improvements are represented by: i) adopted 



state density expression p oc exp (2\/ aU)/U^ (where 
U is the emitting nucleus excitation energy), ii) no 
Maxwellian approximation for energy sampling, 
iii) competition with 7 emission, iv) sub-barrier emis- 
sion. Neutron and proton production are marginally 
affected, while residual nuclei production and alpha 
emission are nicely improved. 

For light residual nuclei, where the excitation en- 
ergy may overwhelm the total binding energy, a sta- 
tistical fragmentation (Fermi Break-up) model is more 
appropriate (see Q, 0, for the FLUKA implemen- 
tation). 

The evaporation/fission/break- up processes repre- 
sent the last stage of a nuclear interaction and are 
responsible for the exact nature of the residuals left 
after the interactions. However, for a coherent self- 
consistent model, the mass spectrum of residuals is 
highly constrained by the excitation energy distribu- 
tion found in the slow stages, which in turn is directly 
related to the amount of primary collisions and follow- 
ing cascading which has taken place in the fast stages. 

IV. NUCLEUS-NUCLEUS COLLISIONS 

The FLUKA implementation of suitable models for 
heavy ion nuclear interactions has reached an opera- 
tional stage. At medium/high energy (above a few 
GeV/n) the DPMJET model is used as described 
in Subsection IIV Al The major task of incorporating 
heavy ion interactions from a few GeV/n down to the 
threshold for inelastic collisions is also progressing and 
promising results have been obtained using a modified 
version of the RQMD-2.4 code (see Subsection llV B|) . 

A. The FLUKA - DPMJET interface 

DPMJET-II.53 [23|, a Monte Carlo model for sam- 
pling h-h, h-A and nucleus-nucleus (A-A) collisions 
at accelerator and cosmic ray energies (E; a fc from 
5-10 GeV/n up to 10 9 -10 n GeV/n) was adapted 
and interfaced to FLUKA. FLUKA implements 
DPMJET-II.53 as an event generator to simulate A- 
A interactions exclusively. DPMJET (as well as the 
FLUKA high energy h-A generator) is based on the 
Dual Parton Model in connection with the Glauber 
formalism. The implementation of DPMJET is also 
considered a possible, future option to extend the 
FLUKA energy limits for hadronic simulations in 
general. 

Internally, DPMJET uses Glauber impact param- 
eter distributions per projectile-target combination. 
These are either computed during initialization of the 
program or can be processed and output in a dedi- 
cated run of DPMJET in advance. The computa- 
tions are CPU intensive for heavier colliding nuclei 
and it would not be practical to produce the required 



MOMT005 



CHEP03, March 24-28, 2003, La Jolla, California 



7 



distributions repeatedly while processing full showers 
in FLUKA. Therefore, a procedure was devised to ef- 
ficiently provide pre-computed impact parameter dis- 
tributions for a complete matrix of projectile-target 
combinations up to a mass number A=246 over the 
whole available energy range |24j . 

FLUKA requires A-A reaction cross sections inter- 
nally in order to select A-A interactions appropriately. 
Hence, a complete matrix of A-A reaction cross sec- 
tions was prepared along with the Glauber impact pa- 
rameter distributions. Owing to the well established 
validity of the Glauber formalism, these cross sections 
can be safely applied down to a projectile kinetic en- 
ergy «1 GeV/n. 

DPMJET is called once per A-A interaction. A 
list of final state particles is returned by DPMJET 
for transport to FLUKA, as well as up to two ex- 
cited residual nuclei with their relevant properties. 
De-excitation and evaporation of the excited residual 
nuclei is performed by calling the FLUKA evapora- 
tion module. 

Work to interface DPMJET-III ;25j is in progress. 
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The FLUKA - RQMD interface 



FIG. 9: The same as fig. for Fe ions. 
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FIG. 8: Double differential neutron yield by 

400 MeV/n Ar ions on thick Al targets. Data are shown 
for six different laboratory emission angles, with the most 
forward on top: histograms are FLUKA results, dots ex- 
perimental data |2fj |. 
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FIG. 10: Fragment charge cross sections for 1.05 GeV/n 
Fe ions on Cu. Stars FLUKA, circles [27l . squares 
(1.5 GeV/n) 0, triangles (1.88 GeV/n) 



Quantum Molecular Dynamics (QMD) approaches 
are a viable solution for A-A reactions. They repre- 
sent an improvement over classical INC codes, thanks 
to their dynamic modelling of the nuclear field among 
nuclcons during the reaction. The treatment of in- 
dividual two-body scattering/interactions is usually 
based on similar approaches for INC and QMD codes. 
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ons. Finally, the remaining energy-momentum con- 
servation issues were resolved taking into account ex- 
perimental binding energies, as in all other FLUKA 
models. After these improvements a meaningful ex- 
citation energy could be computed and the FLUKA 
evaporation model (see Subsection 1111 1J|) was used to 
produce the particles of low energy (in their respec- 
tive rest frames) emitted by the excited projectile and 
target residues. 

Examples of the performances of FLUKA when run 
with the modified RQMD-2.4 code are presented in 
figs. 18191 1101 and ^2 an d compared with experimen- 
tal data, for double differential neutron production at 
400 MeV/n and fragment production at 1.05 GeV/n. 
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PERSPECTIVES 



FIG. 
Al. 



11: The same as fig. [23 for 1.05 GeV/n Fe ions on 



Unfortunately, initialization of the projectile and tar- 
get nuclear states is often difficult and their relativistic 
extension somewhat problematic. 

The RQMD-2.4 [1D,E3 is a relativistic QMD model 
which has been applied successfully to relativistic A- 
A particle production over a wide energy range, from 
« 0.1 GeV/n up to several hundreds of GeV/n. The 
high energy A-A part in FLUKA is already success- 
fully covered by DPMJET. For energies below a few 
GeV/n several models are under development, either 
new (see Section 0), or based on the extension of the 
present intermediate energy h-A model of FLUKA. 
However, a RQMD-2.4 interface has been developed 
meanwhile to enable FLUKA to treat ion interactions 
from sslOO MeV/n up to 5 GeV/n where DPMJET 
starts to be applicable. 

Several important issues had to be addressed. 
RQMD does not identify nuclei in the final state. 
Hence, no low energy de-excitation (evaporation, frag- 
mentation. . . ) is performed for neither the excited 
projectile nor the target residues. This is highly prob- 
lematic, particularly for the projectile residue, since 
its de-excitation usually gives rise to the highest en- 
ergy particle production in the laboratory frame. Se- 
rious energy non-conservation issues are also affecting 
the code, particularly when run in full QMD mode 
(RQMD can be run both in full QMD mode and in 
the so called "fast cascade" mode where it behaves as 
an INC code). Therefore a meaningful calculation of 
residual excitation energies was impossible with orig- 
inal code. 

The adopted solution was to modify the code, re- 
working the nuclear final state out of the available 
information on spectators, correlating the excitation 
energy to the actual depth of holes left by hit nuclc- 
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FIG. 12: 

20 Ne, 40 Ca, 
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Root mean square radii (fm) of selected 
90 Zr initial nuclear configurations versus time 



t = olm/c 

t = 40 fm/c 

t = 80 fm/c 

! = 1 20 fm/c 

I = 1 60 fm/c 

t = 200 fm/c 




FIG. 13: Radial spatial density profile of the same 20 Ne 
configuration shown in fie. 1121 at the beginning of the time 
evolution (solid line) and at subsequent times (At = 40 
fm/c). 
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A new QMD code is also being developed from 
scratch to describe A-A collisions taking into account 
both the effect of stochasting scattering between all 
the nucleons involved and the effect of the nuclear 
potential which acts on each nucleon. Protons and 
neutrons are described as gaussian wave packets of 
fixed widths; the total nuclear wave function is ap- 
proximated by the product of single nucleon wave 
functions. Regarding the Hamiltonian, each group of 
QMD model developers makes its particular choices 
(see, for instance, refs. |M HI IM HI El H3); our 
starting point is a non-relativistic phenomenological 
potential, based on Skyrme interaction, supplemented 
by surface and symmetry terms; we add also the elec- 
tromagnetic repulsion between protons, crucial to de- 
termine low-energy nuclear trajectories. In principle 
one can build an Hamiltonian as sophisticated as de- 
sired to improve the nuclear physical description; in 
practice however one has to meet CPU time require- 
ments. 

The parameters of the model are fixed in order to 
reproduce as accurately as possible the observed nu- 
clear ground state properties. We emphasize that this 
result can be achieved only approximatively, because 
one has to describe with only a few parameters a great 
variety of nuclei, ranging from the lightest to the heav- 
iest ones. For this purpose a gaussian width increasing 
with increasing mass number turns out to be useful, 
as suggested for istance in 35] . On the other hand, 
taking a gaussian width which differs from nucleon to 
nucleon and evolves in time, as done for istance in 
|37|. is practicable only for light nuclei; for the inter- 
mediate mass and heavy ones the required CPU time 
becomes higher and higher. 

We underline that QMD codes are based on Monte 
Carlo simulations: the cross-sections for A-A inter- 
actions are obtained as mean values from hundreds 
and hundreds of events, each of which should involve 
different starting nuclear configurations. In practice, 
it turns out that one can simulate a wide variety of 
different events with only a few initial configurations, 
randomly rotated. 

Reasonable initial configurations can be chosen and 
stored. The time evolution of root mean square radii 



for selected 20 Ne, Ca, 90 Zr configurations is shown 
in fig. El the rms radii slightly oscillate with time be- 
cause the nucleons are not frozen inside the nucleus: 
each of them moves in the potential well originated 
from all the others because of its Fermi momentum, 
different from zero. Note that QMD initial configu- 
rations are not classical Hamiltonian minima, which 
would break the Pauli principle Q • One has to make 
sure that no nucleons escape at least for a time of the 
order of that required to A-A collisions to take place 
(~ hundreds of fm/c); only configurations which do 
not originate this spurious emission are selected and 
stored for subsequent simulations. . 

The stored configurations should also provide rea- 
sonable values for density and momentum distribu- 
tions. 

The radial spatial density profile is plotted in fig. 1131 
for the same 20 Ne nucleus whose root mean square 
radius evolution is shown in fig. 1121 each curve refers 
to a different time during evolution. One can see that, 
as well as the radius, also the density oscillates around 
its typical mean value. Similar plots can be obtained 
for the time evolution of the momentum distribution. 

It is planned to couple the dynamical nuclear evo- 
lution predicted by our QMD, which gives the de- 
scription of the first stage of the reaction, with the 
FLUKA preequilibrium module, to describe the de- 
excitation of the fragments formed and to study 
deeply the fragmentation process, implementing suit- 
able models. 

Moreover, a promising task is represented by the 
coupling of FLUKA with a Monte Carlo code [33] de- 
veloped at Milan University and based on Boltzmann 
Master Equation theory, as a tool to treat ion-ion in- 
teractions below 100 MeV/n. 
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